Semiclassical features of rotational ground bands 

A. A. Raduta"' 6 ), R. Budaca a ^ and Amand Faessler c ) 
^Institute of Physics and Nuclear Engineering, Bucharest, POB MG6, Romania 
b > Academy of Romanian Scientists, 54 Splaiul Independents, Bucharest 050094, Romania and 
c ) Institut filr Theoretische Physik der Universitat Tubingen, 
q ; Auf der Morgenstelle 14, D-72076 Tubingen, Germany 

(Dated: April 8, 2010) 

; Abstract 

A time dependent variational principle is used to dequantize a second order quadrupole bo- 
son Hamiltonian. The classical equations for the generalized coordinate and the constraint for 
angular momentum are quantized and then analytically solved. A generalized Holmberg-Lipas 
formula for energies is obtained. A similar J (J + 1) dependence is provided by the coherent 
state model (CSM) in the large deformation regime, by using an expansion in powers of 1/x for 



O 
=5 



energies, with x denoting a deformation parameter squared. A simple compact expression is also 



o 



possible for the near vibrational regime. These three expressions have been used for 44 nuclei 
covering regions characterized by different dynamic symmetries or in other words belonging to 
the all known nuclear phases. Nuclei satisfying the specific symmetries of the critical point in 



o 



X 



the phase transitions 0(6) — > SU(3>), SU(5) — > SU (3) have been also considered. The agreement 
between the results and the corresponding experimental data is very good. This is reflected in 
very small r.m.s. values of deviations. 
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I. INTRODUCTION 



One of the big merits of the liquid drop model consists of that it defines in a consistent 
way the rotational bands. Many theoretical efforts have been made for the description of 
excitation energies as well as of electromagnetic transitions probabilities. One of aims, 
in the beginning era, was to obtain a closed formula for the ground band energies which 
explain the deviations from the J( J + 1) pattern. Various methods have been proposed 
which were mainly based on the variational moment of inertia principle [lHSfl. These 
approaches proposed for ground band energies a series expansion in terms of J( J + 1) 
term. The weak point of these expansions is that they do not converge for high angular 
momenta. The first attempt to avoid this difficulty was due Holmberg and Lipas [4] who 
proposed a square root of a linear expression of J( J+ 1). This expression proves to work 
better than a quadratic expression in J(J + 1). 

Here we address the question whether this formula can be improved such that it 
could be extended to the region of states with high angular momenta. In the present 
paper we offer three solutions for this problem, each of them being obtained in a distinct 
manner. One solution is based on a semiclassical treatment of a second order quadrupole 
boson Hamiltonian. The remaining two expressions for the ground band energies are 
given by asymptotic and near vibrational expansions respectively, of an angular mo- 
mentum projection formula. The three expressions obtained for energies are used for a 
large number of nuclei. The above sketched project has been achieved according to the 
following plan. In Section II we present the semiclassical approach in connection with a 
quadratic quadrupole boson Hamiltonian. In Section III the angular momentum projec- 
tion method is described. Numerical applications are presented in Section IV, while the 
final conclusions are drawn in Section V. 
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II. SEMICLASSICAL TREATMENT OF A SECOND ORDER QUADRUPOLE BOSON 
HAMILTONIAN 



For a moment we consider the simplest quadrupole boson (b\ , — 2 < \i < 2) Hamil- 
tonian: 

h = a x x; ^ + a 2 j2( b l b -» + M^x-r- ci) 

Here we are interested to study the classical equations provided by the time dependent 
variational principle associated to H: 

5 J (V\H - ih^\9)dt = 0. (2.2) 

If the variational states span the whole Hilbert space of boson states, then solving the 
variational equations is equivalent to solving the time dependent Schrodinger equation 
which is in general a difficult task. Therefore we restrict the trial function to a coherent 
state which we hope is the suitable state for describing the semiclassical feature of the 
chosen system: 

|^) = Exp \z Q bl - z* b + z 2 (b\ + bl 2 ) - z* 2 (b 2 + 6-2)] • (2.3) 
Indeed the coherence property results from the obvious equation satisfied by 

bM = (<Wo + (^2 + 5^ 2 )z 2 )\tP). (2.4) 



In order to write explicitly the equations emerging from (|2.2[) we have to calculate first 
the averages of H 

H=(^\H\9), (2.5) 

as well as of the action operator —ifr-ii- The variational equation (|2.2[) yields the following 
classical equations for the complex coordinates z k and z* k : 

du . t ., an . t . 
dz- = -' h ^ 84= °' 
dH dH 

— = -2ihz 2 , — = 2ihz 2 . (2.6) 
oz 2 oz\ 

Note that the coordinates z k and z* k define a classical phase space while H plays the role 
of a classical Hamilton function. For what follows it is useful to bring these equation to 



a canonical form. This is achieved by the transformation: 

g _ 2 ( fc+2) /4 Re( ^ ); pi = h2 (^2)/ Hm{zkl fc = 0)2? i = h±l. (2.7) 

Indeed, in the new coordinates the classical equations of motion become: 

&H . &H . 

-w- = -Pk, tt~ = Qk- (2.8) 

oq k op k 

In terms of the new coordinates, the Hamilton function is written as: 

A l + 2A 2 ( 2|2 ,, A l -2A 2 ( 2 2 ^ 
^ = ^ Wi + feH — ^ — (P1+P2) 

= y(rf + d) + ^(p?+rf), (2-9) 

where we denoted by A = A x + 2A 2 and A' = A\ — 2A 2 . Eqs. (|2.8)> provide the connection 
between the generalized momenta and the coordinates time derivatives: 

Taking into account these relations, the classical energy function becomes: 

n = ^(ql + q 2 2 ) + ^( q 2 + q l). (2.11) 

For what follows it is useful to use the polar coordinates: 

qi=rcos9, q 2 = r sin 9, (2.12) 

for the Hamilton function: 

H = ^ 7 (r 2 + r 2 ^) + |r 2 . (2.13) 

The classical system described by H is exactly solvable since the number of degrees of 
freedom is equal to the number of constants of motion. Indeed, taking the time deriva- 
tives of % and £3, the third component of a pseudo-angular momentum acting in a 
fictitious boson space, one obtains: 

H = 0. £3 = 0. (2.14) 

The components of the pseudo-angular momentum are defined in Appendix A. Here 
we need the conserved component: 

£3 = ^{<HP2 -92Pi)- (2.15) 



Its constant value is conventionally taken to be: 



(2.16) 



which allows us to express the angular variable derivative in terms of the radial one: 
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2A'L 



fir 



2 • 



Thus, the energy function written in the reduced space, becomes: 



U 



h 2 



:r 2 + 



2A'L 2 A 



h 2 



2A r ' r 2 2 2A 
We recognize in the effective potential energy: 

2A'L 2 A 



-r 2 + V eff (r). 



Veff(r) 



+ 77 r » 



(2.17) 



(2.18) 



(2.19) 



r* 2 
just the Davidson potential. 

Instead of solving the classical trajectories and then quantizing them, here we first 
quantize the energy by replacing 



h 2 r d 



Thus, one arrives at the Schrodinger equation: 

A' d 2 2A'L 2 A 



+ 



_l r ^ 

2 2 



2 dr 2 r 

Making use of the change of variable and function: 

^4 



u(r) = eu(r). 



(2.20) 



(2.21) 



x 




-r , u(r) = e 2x s f(x), 



one obtains the following differential equation: 



d 2 



x 



dx 2 + [ 2S+ 2- X )d-x" 



d 



2s 2 - s - 2L 2 



+ 



— s 



(2.22) 



f(x) = 0. (2.23) 



2x ' 2y/AA 4 

This should be compared with the differential equation for the Laguerre polynomials: 



d 2 , J 
x ^r^ + {m + 1 - x)— + n 
ox 1 ox 



L™(x) = 0. 



(2.24) 



Indeed, the two equations are identical provided the following equations hold: 



From the last equation we derive the expression of s as a function of L. The positive 
solution is: 

s = ^(l + v / l + 16L 2 ). (2.26) 
The second equation (|2.25|) yields for the energy e the following expression: 



e = 2sJ {A\ - 4A|) (n + X - + -VT+ 16L^j , n = 0, 1, 2, ... L = 0, 1, 2, ... (2.27) 

An approximative expression may be obtained by expanding first the Davidson po- 
tential V e ff around its minimum r given by the equation: 

rl = 2L^, (2.28) 

and truncating the expansion at the quadratic term. The result for the energy function 
is: 

h 2 

H = —r 2 + 2A(r - r ) 2 + 2Lv / Z4 7 . (2.29) 

Quantizing this Hamilton function we obtain an eigenvalue equation for a harmonic 
oscillator whose energy is: 

E nL = 2v / A4 7 (n + + 2Ly/AA' = 2^(A 2 -Aa 2 ) (n + ^ + Lj , n = 0, 1, 2, ... (2.30) 

We remark the fact that the two spectra coincide when L is large: 

E n L ~ en,L, for L = large. (2.31) 

Note that the initial boson Hamiltonian could be easily diagonalized by a suitable 
chosen canonical transformation: 

K = ubi - vi-rb-v 

K = Ub, - V(-ybi„. (2.32) 
Indeed, the coefficients U and V may be chosen such that: 



MI' 



E%. (2.33) 
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The second equation provides a homogeneous system of equations for the tranformation 
coefficients 

/ A 1 2A 2 \ f u\ ( u\ 

\=E\ , (2.34) 
-2A 2 -A l J \V J \V J 

which determine U and V up to a multiplicative constant which is fixed by the first 
equation which gives: 

U 2 - V 2 = 1. (2.35) 



The compatibility condition for Eq. (|2.34[) gives E = \A4f — 4 A\, and therefore the eigen- 
values of H are: 



E n = ^A 2 -AA 2 [n + -j. (2.36) 

The frequency obtained is half the one obtained through the semiclassical approach. The 
reason is that here the frequency is associated to each of the 5 degrees of freedom while 
semiclassically the frequency is characterizing a plane oscillator. Note that the pseudo- 
angular momentum L is different from the angular momentum in the laboratory frame 
describing rotations in the quadrupole boson space: 



The expected value of the angular momentum square is 



J M = V6 (ftj&a) . (2.37) 



V\J 2 \V) = 2 



ii + il + ^ (pi + pI) 



(2.38) 



Since the variational function |^) is not eigenstate of J 2 , the above mentioned average 
value is not a constant of motion. Indeed, it is easy to check that: 

[ g t l ; = ^ (A '- A) ( gm + q 2P2 ) ? 0. (2.39) 

It is instructive to see whether we could crank the system so that the magnitude of an- 
gular momentum is preserved, i.e. 

(^|J 2 |^) = h 2 J(J+ 1). (2.40) 

Using the polar coordinates the above equation becomes: 

3 -Lf> + 1 ^- + ^ = J(J+l). (2.41) 



This equation is treated similarly with the energy equation. Thus by the quantization: 

A' 2 8 
at or 



Eq. (|2.41|) becomes a differential equation for the wave function describing the angular 
momentum: 

d 2 $ / 4L 2 r 2 \ + 

+ ( + — ) $ = \ — (2.43) 



<9r 2 V yl' 2 r 2 A' 2 I 
Making the change of variable and function: 

x = — , $ = e"^x s \l/, 

we obtain the following equation for \1/: 



3A' 2 



(2.44) 



<9 2 # / 1\ d^> 

dx 2 \ 2 / dx 



2s 2 -s- 2L " 



2x 



12A' 



s - - 1 * = 0. 



(2.45) 



This equation admits the Laguerre polynomials (x) with the quantum numbers de- 
termined as follows: 



1 11/ 16L 2 JJ+1 

m =2s , s = - + -a/H 9-, — - — rr^ 

2' 4 4 V A' 12A' 



, 1 1 / 16L 2 

n H h -a/1 H o-- 

2 4 V A' 



(2.46) 



The last relation (|2.46|) can be viewed as an equation determining L: 



L 



12 v 2' 



n i/2 



(2.47) 



On the other hand taking the harmonic approximation for the potential term in Eq. (|2.41[) 
one obtains the classical equation for a harmonic oscillator from which we get: 

J(J + 1) = 12A' (n' + + 12L. (2.48) 

Reversing this equation one can express the pseudo-angular momentum L in terms of 
the angular momentum J: 



= 12 [ n+ 2 ] 



(2.49) 



Replacing, successively the expressions for L, (|2.47H2.49|) , into energy equations (|2.27l 
|2.30|) , we obtain four distinct expressions for the energies characterizing the starting 
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Hamiltonian H. 



E^J = ^A 



E^j = VX4 7 



2,/ + l + J(J . + 1) -A / (2n / +l) 



6 



2n + 1 + 2\ 



A,{ , 1 



E^j = VAA 



2n + 1 + -\- I 



6 



A' (2ri + 1) 



E ( 3j = ^AA 



2n + l 



1 + 4 



J(J+l) 
6 



1 2 



A' (2n' + r 



-(A') 



(2.50) 
(2.51) 

(2.52) 

(2.53) 



Remark the fact that for a fixed pair of (n, n') each of the above equations define a ro- 
tational band: The lowest band corresponds to (n, n') = (0,0) and defines the ground 
band. Except for the band energies E^j which exhibits a J(J+ 1) pattern the other three 
bands have the same generic expressions. Thus the excitation energies have the form: 



Ej = a 



^/l + bJ(J + l) + cJ 2 (J+l) 2 -l 



(2.54) 



which is a generalization of the Holmberg-Lipas formula ||4j]. 

We recall that we required that the average value of J 2 equals h 2 J(J + 1). Subse- 
quently we eliminated the energy dependence on the pseudo-angular momentum L. In 
this way we projected approximately the angular momentum from the variational state. 
In the next section we shall show that the exact treatment of the angular momentum 
projection yields also a closed formula for energy as function of J( J + 1). 



III. THE METHOD OF ANGULAR MOMENTUM PROJECTED STATE 



For the sake of simplicity here we consider a simple form for the variational state 



l**> = e 



<2(&L-6 2 o) 



10} 



(3.1) 



in connection with the following quadrupole boson Hamiltonian: 



H = A 1 J2blh fl + A 2 J 2 



(3.2) 
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The vacumm state for the quadrupole boson operators is denoted by |0) while d is a real 
quantity which plays the role of the deformation parameter. The reason is the fact that 
the average value of the quadrupole moment, written in the lowest order in terms of 
quadrupole boson operator, with the function is proportional to d. The component 
of a given angular momentum is obtained by a projection procedure: 

¥>Kr = nPp&Pi, (3-3) 
where Pm K denotes the angular momentum projection operator: 

p mk = J Djz K (n)R(n)dn, (3.4) 

with Dj£ K denoting the Wigner functions and R(Q) a rotation defined by the Euler angle 
f2. The system energy is defined as the average value of H with the projected state: 



J?(d^ 
lf\d* 



Ef ee (vf M \HW%) = Md 2 ^^ + A 2 J(J + 1), (3,5) 

where we denoted by if^ the overlap integral: 

lf\x)=2 [ Pj(y)e xP2iy) dy, x = d 2 . (3.6) 
Jo 

The k derivative of this integral is denoted by: 

tf'M - & (3-7) 

The normalization constant for the projected state has the expression: 

i^Nfy 2 = {2J + l)lfe- d2 . (3.8) 

These integrals have been analytically calculated in Ref.H- Actually the energies pre- 
sented here refer to the ground band described by the coherent state model (CSM) which 
considers simultaneously three interacting bands, ground, beta and gamma. In the 
asymptotic limit of the deformation parameter d the ground band energies have the 
expression IS] 



E 



(g,asym) A\ 



2 



+ A 2 J(J + 1), (3.9) 
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with 



G, 



9 



x(x — 2) + ( J 



+ 



3.r 



10 13 
1 + — + 



3x 



x- 



4 

_ 9x 
J(J+1) 



10 37 

3+ — + 



X 



X, 



9x 3 



J 2 (J + 1) 



X 



d 2 . 



(3.10) 



It is worth to mention that Eq. (|3.9|) is similar to the generalized HL formula, with the 
difference that here the coefficients of the terms J( J + 1) and J 2 (J + l) 2 have explicit 
expressions in x. Moreover there appears an additional J( J + 1) term outside the square 
root symbol. The expression (|3.9|) is obtained by replacing the series expansion in 1/x, 
associated to the ratio x^ 



(i) 
j 

(0)/ 



x- 



r(D 
■7 

r(0) 



1 

1 

3x 



5 
9X 2 



37 
27x3 



1 5 
6x 18x 2 



+ 



13 

18x 3 



J(J+l) 



54a; 3 



f(J+lf + 0(x- 



(3.11) 



by a faster convergent one. 

According to Ref-Ba], for the near vibrational regime (d- close to zero) the ground state 
band energies have the expressions: 



E 9 j 



vib 



+ 



A 1 
27 



J 



J 



9 (J+l)(J + 2) 

2 ^ 2(2J + 3) X ^ 2 (2J + 3) 2 (2J + 5) 
(J+l)(J + 2) ^ 



x 



-X 



+ A 2 J(J + 1). 



(3.12) 



2 (2J + 3) 3 (2J + 5)(2J + 7)' 

For the sake of completeness we present the derivation of the two expressions for the 
ground band energies, in the rotational and near vibrational limits, in Appendix B. 



IV. NUMERICAL RESULTS 



Since the expressions (|3.9[), (|3.11|) and (|3.12|) are based on series expansion in 1/x and 



x, respectively, it is worth showing how far are the truncated expansions from the exact 

7 (i) 

energies. Aiming at this goal in Fig. 1 and Fig. 2 we plotted the ratio d -fa and the asso- 
ciated truncated series for large and small values of d respectively, as functions of d for 
two angular momenta: J = 12 and J = 16. In the case of asymptotic regime we consid- 
ered also the square root expression. In this case one defines an existence interval of d for 
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which Gj > 0. The lower bounds of these intervals for J running from to 30 are listed 
in Table I. From Fig. 1 we see that for d > 3 the used expressions for energies achieve 
the convergence even for high angular momenta. Concerning the energies for the near 
vibrational regime one notes that we use a power series of x and therefore one may think 
that such an expansion is valid for x < 1. However, we notice that the coefficients of this 
expansion are depending on J and moreover are under unity numbers. The larger J the 
smaller are these coefficients. This fact infers that the convergence radius is larger than 
unity and is an increasing function of the angular momentum. As a matter of fact this 
is confirmed in the plot shown in Fig. 2. Comparing the curves from Figs. 1 and 2 one 
may say that there is a small interval of d were the asymptotic and small x expansions 
are matched. This allows us to assert that the reunion of the two formulas, (|3.9|) and 
(|3.12|) , assures an overall description of nuclei ranging from small to large deformation. 
In Fig. 3 we plotted the term Gj involved in the energy expression (|3.9|) as a function of 
the deformation parameter d. Except for J = and J = 2 all the other functions vanish 
for a specific value of d which are, in fact, the lower bounds of the existence interval. 



The basic expressions for energies (|2.54|) , (|3.9|) and (|3.12|) have been used for a large 
number of nuclei grouped according to the nuclear phase to which they belong. Thus, 
for well deformed nuclei behaving like axially deformed rotator the ratio E 4 +/E 2 + 
should be close to the value of 3.3 while for the near vibrational region one expects a 
ratio close to the value 2. Between these two extreme values are placed gamma unstable 
nuclei where the ratio may run in the interval of 2.5-3.0. The deviation from axial sym- 
metry can affect the ratio mentioned above. Thus 228 Th exhibits some specific feature of a 
triaxial nucleus with an equilibrium value 7 = 30°. The corresponding ratio E 4 +J E 2 + is 
equal to 3.24. According to the IBA (Interacting Boson Approximation) model [Zr 3l the 
nuclei belonging to the three groups mentioned above are described by the irreducible 
representations of some dynamic groups as SU (3), SU(5) and 0(6). Since the nuclei de- 
scribed by a certain symmetry group exhibit some specific distinct properties one says 
that these form a certain nuclear phase. According to Casten Igfl all nuclei of the periodic 
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25 



20 



15 
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• Exact 

Truncated series 
Square root expression 





FIG. 1: d 2 Ij /I j is plotted as a function of d for two values of angular momentum. Two 
approximations of this function are also presented. One is a truncated expansion in 1/x, while 
the other one is given by a square root expression which is slightly faster convergent than the 
previously mentioned expansion. 13 





FIG. 2: d 2 Ij /Ij is plotted as a function of d for two values of angular momentum. This is 
compared with the function given by the near vibrational approximation from Eq. (|3.i2[) . 
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2 

d 



FIG. 3: Gj, is plotted as function of d for some angular momenta. Note that excepting the cases 
of J = 2, 4 all other functions Gj get negative for d smaller than a critical value. These limiting 
values are listed in Table I. 

table may be placed on the sides of a triangle having in vertexes the three symmetries 
mentioned above. On each side which links two adjacent symmetries one expects a crit- 



ical transition point between the two adjacent phases. Few years ago, Iachello [10, [Hi] 
advanced the idea that each of the critical nuclei laying on the three triangle sides cor- 
respond to specific symmetries. Thus, the transition 0(6) — > SU(5) is characterized by 
a critical symmetry which is E(5). Representatives for the E(h) symmetry are 104 Ru and 
102 Pd characterized by specific ratios E 4 +/E 2 + = 2.48, 2.29, respectively. In the transition 
SU(5) SU(3), the critical point is close to 3. Such nuclei are 150 Nd, 152 Sm, 154 Gd and 
156 Dy. Indeed, they prove to be critical points for the mentioned phase transition when 
the entire isotopic chains are considered. 

The question to be answered is whether the compact energy formulas obtained in this 
paper are able to describe the ground band energies for all nuclei mentioned above. The 



theoretical values for energies labeled by Th(l) are obtained with Eq. (|3.9|) if d is large 
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or with Eq. (|3.12|) when d is smaller than 2. The calculated energies labeled by Th(2) are 
obtained with Eq. (|2.54[) . The parameters A\, A 2 , d were obtained by a least mean square 
fitting procedure while a, b, c by fixing the energies of three particular levels. 

The agreement of calculated and experimental excitation energies is judged by the 
r.m.s (root mean square) values of the deviation, denoted by 



X 



\ 



N {Efh _ E E XP)2 

2^ n 



(4.1) 



The fitting procedure yields for the coefficients b and c double precision numbers, which 
are presented, in tables, in a truncated form. Since the square root formula provides 
energies which are quite sensitive to small variations for the parameters b and c, we give 
their values with a suitable large number of digits. Indeed, with the listed parameters we 
get the energies corresponding to the exact parameters yielded by the fitting procedure. 
Comparing the values of c for different nuclei, one remarks that the parameter acquires 
larger values for smaller deformation parameter. For two nuclei, 248 Cm and 180 Os, the 
parameter c gets negative values which annihilates a part of the contribution coming 
from the J (J + 1) term. 



J 



2 4 6 



10 12 14 16 18 20 22 24 26 28 30 



1.55 0.65 1.21 1.43 1.59 1.71 1.82 1.91 1.99 2.07 2.14 2.21 2.27 2.33 



TABLE I: The smallest value of d, for which Gj is positive. 
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22S Th 

Exp. Th(l) Th(2) 



Exp. 



>,32 Th 

Th(l) Th(2) 



Exp. 



!32 v 

Th(l) Th(2) 



Exp. 



Th(l) Th(2) 



236 ;j 

Exp. Th(l) Th(2) 



Exp. 



:38 r/ 

Th(l) Th(2) 



2+ 
4+ 
6+ 
8+ 
10+ 
12+ 
14+ 
16+ 
18+ 
20+ 
22+ 
24+ 
26+ 
28+ 
30+ 



57.76 

186.82 

378.18 

622.50 

911.80 

1239.4 

1599.5 

1988.1 

2407.9 



57.61 

186.90 

378.51 

622.73 

911.49 

1238.58 

1599.27 

1989.81 

2407.05 

2848.18 

3310.53 

3791.4 

4287.89 

4796.81 

5314.41 



57.485 

186.82 

378.78 

623.29 

911.80 

1237.99 

1597.63 

1988.1 

2407.9 

2856.32 

3333.13 

3838.43 

4372.51 

4935.76 

5528.64 



49.37 

162.12 

333.2 

556.90 

827.0 

1137.1 

1482.8 

1858.6 

2262.9 

2691.5 

3144.2 

3619.6 

4116.2 

(4631.8) 

(5162) 



49.66 

162.966 

334.64 

558.25 

827.44 

1136.62 

1481.12 

1857.14 

2261.56 

2691.82 

3145.71 

3621.32 

4116.88 

4630.77 

5161.41 



49.31 

162.12 

333.63 

557.65 

827.69 

1137.77 

1482.77 

1858.6 

2262.1 

2690.9 

3143.29 

3618.06 

4114.4 

4631.8 

5169.96 



47.57 

156.57 

322.60 

541.00 

805.80 

1111.5 

1453.7 

1828.1 

(2231.5) 

(2659.7) 



47.46 

156.34 

322.59 

540.94 

805.75 

1111.64 

1453.81 

1828.13 

2231.13 

2659.9 

3111.96 

3585.21 

4077.8 

4588.11 

5114.65 



47.54 

156.57 

322.96 

541.33 

805.97 

1111.5 

1453.3 

1827.61 

2231.5 

2662.76 

3119.8 

3601.5 

4107.12 

4636.2 

5188.48 



43.50 

143.35 

296.07 

497.04 

741.2 

1023.8 

1340.8 

1687.8 

2063.0 

2464.2 

2889.7 

3339 

3808 

(4297) 



43.72 

143.93 

296.84 

497.63 

741.3 

1023.38 

1339.79 

1687.24 

2062.99 

2464.74 

2890.57 

3338.8 

3807.96 

4296.74 

4803.88 



43.49 

143.35 

296.07 

497.02 

741.15 

1023.69 

1340.41 

1687.8 

2063.08 

2464.12 

2889.32 

3337.54 

3808.0 

4300.17 

4813.74 



45.24 
149.48 
309.78 
522.24 
782.3 
1085.3 
1426.3 
1800.9 
2203.9 
2631.7 
3081.2 
(3550) 
(4039) 
(4549) 
(5077) 



45.516 

150.33 

311.34 

524.29 

784.32 

1086.49 

1426.07 

1798.79 

2200.85 

2628.97 

3080.31 

3552.43 

4043.2 

4550.77 

5073.52 



45.22 

149.48 

309.97 

522.724 

783.07 

1086.07 

1426.86 

1800.9 

2204.11 

2632.91 

3084.22 

3555.4 

4044.28 

4549.0 

5068.04 



44.92 

148.38 

307.18 

518.10 

775.9 

1076.7 

1415.5 

1788.4 

2191.1 

2619.1 

3068.1 

3535.3 

4018.1 

4517 

5035 



44.73 

147.93 

306.90 

517.8 

776.26 

1077.35 

1416.33 

1788.68 

2190.24 

2617.29 

3066.53 

3535.05 

4020.31 

4520.08 

5032.4 



44.88 

148.38 

307.71 

518.96 

777.492 

1078.41 

1416.88 

1788.4 

2188.89 

2614.76 

3062.93 

3530.78 

4016.08 

4517.0 

5031.98 



E 



3.23 



3.28 



3.29 



3.30 



3.30 



3.30 



0.56 



0.6 



1.13 



2.2 



0.17 



1.00 



0.53 



0.95 



2.21 



3.16 



1.40 



2.48 



182.8720 

4.0861 

2.7062 



1909.4426 
0.0101534 
5.4707523-10" 6 



233.4401 

3.2430 

3.0822 



3339.3037 
0.00495222 
1.026518-10- 6 



298.6610 

2.6094 

3.3385 



3472.2796 
0.00458875 
1.102772-10" 6 



221.7839 

2.9879 

3.2192 



3443.7877 
0.00423038 
1.425528- 10~ 7 



386.4548 

1.9129 

3.6189 



5932.3853 
0.00254985 
7.754542 -lO" 8 



502.4456 

1.0667 

3.8524 



5945.5847 
0.00252536 
6.979548- 10" 8 



TABLE.,: Expen.en.al (Ex., and theoretical M and Th(2)> excitation enemies, for several nuclei, -Th llZD, -Th M, «U Uffl, 
234 U (ISO, 236 U (20D, 238 U 121D, are given in units of keV. The predictions labeled by Th(l) are obtained with the square root formula 
given by the asymptotic expansion of the CSM-ground band energies while Th(2) are obtained by the generalized HL expression. The 
parameters for Th(l) calculations, i.e. A\, A2, d were obtained by a least mean square procedure while those of the set Th(2), a, b, c, by 
fixing three particular energy levels, which are underlined. The obtained parameters are also listed. To have a hint about the agreement 
between the theoretical and experimental excitation energies, for each case the r.m.s value of discrepancies, denoted by x> is a l so given. 
The values of A\, A2, a, x are given in keV while, d, b, c are dimensionless. Having in view a possible classification of the considered 
nuclei, the ratio /Eg + is also given. 
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236p u 


238 Pu 


24 °Pu 


242 Pu 


248 Cm 


J* 


Exp. 


Th(l) 


Th(2) 


Exp. 


Th(l) 


Th(2) 


Exp. 


Th(l) 


Th(2) 


Exp. 


Th(l) 


Th(2) 


Exp. 


Th(l) 


Th(2) 


2+ 


44.63 


44.514 


44.602 


44.076 


43.818 


44.06 


42.824 


42.88 


42.82 


44.54 


44.14 


44.49 


43.40 


43.41 


43.33 


4+ 


147.45 


147.23 


147.45 


145.95 


145.25 


145.95 


141.69 


141.84 


141.69 


147.3 


146.29 


147.3 


143.6 


143.80 


145.6 


6+ 


305.80 


305.55 


305.82 


303.38 


302.44 


303.61 


294.32 


294.52 


249.34 


306.4 


304.52 


306.20 


298.1 


299.18 


298.88 


8+ 


515.7 


515.80 


515.939 


513.58 


512.76 


514.16 


497.52 


497.73 


497.61 


518.1 


516.00 


517.98 


505.0 


506.59 


506.30 


10 


+ 


773.5 


773.6 


773.5 


773.48 


773.027 


774.231 


(747.8) 


747.93 


747.93 


778.6 


777.19 


778.86 


760.7 


762.36 


762.29 


12 


+ 


1074.3 


1074.31 


1074.1 


1080.1 


1079.84 


1080.33 


(1041.8) 


1041.61 


1041.71 


1084.4 


1084.14 


1084.78 


1061.3 


1062.39 


1062.77 


14 


+ 


1413.6 


1413.31 


1413.6 


1429.1 


1429.73 


1429.1 


(1375.6) 


1375.54 


1375.6 


1431.7 


1432.67 


1431.7 


1402.5 


1402.35 


1403.35 


16 


+ 


1786.0 


1786.21 


1788.29 


1818.5 


1819.41 


1817.49 


(1746.9) 


1746.84 


1746.67 


1816.7 


1818.59 


1815.83 


1779.6 


1777.89 


1779.6 


18 


+ 




2188.99 


2195.01 


2244.9 


2245.79 


2242.86 


(2153.1) 


2153.04 


2152.48 


2236.0 


2237.82 


2233.73 


2187.7 


2184.82 


2187.09 


20+ 




2618.05 


2631.13 


2705.7 


2706.08 


2703.0 


(2591.9) 


2592.01 


2591.07 


2686 


2686.51 


2682.38 


2621.5 


2619.16 


2621.59 


22+ 




3070.18 


3094.55 


3198.8 


3197.78 


3196.12 


(3062.2) 


3061.95 


3060.93 


3163 


3161.07 


3159.2 


3077.2 


3077.22 


3079.06 


24+ 




3552.59 


3583.61 


3720.8 


3718.68 


3720.8 


(3560.9) 


3561.32 


3560.9 


3662 


3658.24 


3662.0 


3552.4 


3555.65 


3555.75 


26+ 




4032.81 


4097.01 


4265.2 


4266.82 


4275.97 


(4089) 


4088.78 


4090.18 


4172 


4175.02 


4188.99 


4048.2 


4051.4 


4048.2 


28 


+ 




4538.68 


4633.81 




4840.48 


4860.82 




4643.16 


4648.2 




4708.76 


4738.69 


4564.5 


4561.71 


4553.19 


30+ 




5058.3 


5193.28 




5438.15 


5474.79 




5223.42 


5234.61 




5257.04 


5309.9 




5084.14 


5067.79 


t + 


3.30 






3.31 






3.31 






3.31 






3.31 






X 






0.18 


1.72 




0.99 


3.24 




0.19 


0.57 




1.89 


4.98 




1.95 


3.28 


Ai 


a 


486.0023 


4857.3300 


397.4155 


5443.1075 


232.3635 


4487.9012 


692.5499 


6375.5552 


705.2236 


10281.7179 


M 


b 


1.2356 


0.00307187 


2.4878 


0.00270563 


3.42535 


0.0031902 


0.3216 


0.00233268 


0.0133 


0.00140828 


d 


c 


3.8514 


4.909159 10 -7 


3.9353 


5.862828-10~ 7 


3.4867 


8.397898-10" 7 


4.2578 


2.196286- 10" 7 


4.2426 


-9.361609-10" 8 



TABLE III: The same as in Table II but for a different set of nuclei: 236 Pu [20], 238 Pu |21fl, 240 Pu |22D, 242 Pu I23Q, 248 Cm 





154 Nd 


156 Nd 


1S6 Sm 


158 Sm 


160 Gd 


162 Gd 




Exp. 


Th(l) 


Th(2) 


Exp. 


Th(l) 


Th(2) 


Exp. 


Th(l) 


Th(2) 


Exp. 


Th(l) 


Th(2) 


Exp. 


Th(l) 


Th(2) 


Exp. 


Th(l) 


Th(2) 


2+ 


70.8 


70.55 


70.70 


66.9 


67.13 


67.06 


75.89 


75.60 


75.633 


(72.8) 


72.75 


72.581 


75.26 


75.153 


75.22 


72.1 


71.73 


71.89 


4+ 


233.2 


233.80 


233.2 


221.8 


221.93 


221.8 


249.71 


249.64 


249.71 


(240.3) 


240.538 


240.3 


248.52 


248.34 


248.52 


237.3 


236.94 


237.3 


6+ 


481.9 


481.65 


282.13 


460.4 


460.48 


460.42 


517.07 


517.13 


517.21 


(498.4) 


499.16 


499.41 


514.75 


514.77 


514.02 


490.8 


490.93 


491.25 


8+ 


810.1 


810.26 


810.48 


777.9 


777.76 


777.9 


871.9 


871.90 


871.9 


(844.5) 


843.26 


844.5 


867.9 


868.0 


868.15 


827.3 


827.32 


827.3 


10+ 


1210.8 


1211.23 


1210.8 


1168.9 


1168.57 


1168.87 


1307.4 


1307.72 


1307.67 


(1266.7) 


1267.24 


1269.02 


1300.7 


1300.84 


1300.7 


1238.9 


1239.09 


1238.75 


12+ 


1677.3 


1677.3 


1676.12 


1628.4 


1628.05 


1628.23 


1819.3 


1818.97 


1819.3 


(1765.8) 


1765.92 


1765.8 


1806.3 


1806.12 


1805.6 


1719.5 


1719.26 


1719.5 


14+ 


2202.4 


2201.89 


2200.46 


2151.6 


2151.99 


2151.6 


2400.8 


2400.87 


2402.81 


(2334.9) 


2334.77 


2327.52 


2377.3 


2377.21 


2376.51 


2261.3 


2261.39 


2264.46 


16+ 


2779.0 


2779.24 


2779.0 


2737.0 


2736.83 


2735.47 




3049.5 


3055.43 




2969.98 


2947.08 


3008.1 


3008.19 


3008.1 




2859.77 


2869.72 


18 


+ 


(3399.3)? 


3404.46 


3408.08 




3379.62 


3377.18 




3761.62 


3775.4 




3668.38 


3617.8 




3693.96 


3696.07 




3509.48 


3532.43 


20+ 




4073.36 


4085.0 




4077.89 


4074.82 




4534.54 


4561.74 




4427.32 


4333.63 




4430.11 


4437.09 




4206.31 


4250.64 


K 
El 


+ 


3.29 






3.32 






3.29 






3.30 






3.30 






3.29 






X 




0.32 


0.82 




0.23 


0.57 




0.21 


0.78 




0.60 


2.95 




0.12 


0.40 




0.23 


1.21 


A x 


a 


441.8371 


5240.0423 


241.6860 


5338.3543 


251.2029 


4267.15 


259.4635 


12111.9787 


485.7106 


6680.1981 


460.9196 


4947.3299 


A 2 


b 


4.3446 


0.00451888 


6.5982 


0.00420209 


7.4572 


0.00593371 


7.2031 


0.00200342 


4.8487 


0.00376825 


4.5633 


0.00486469 


d 


c 


3.4184 


1.524693-10~ 6 


3.24587 


1.952257- 10" 6 


3.1471 


4.471374-10 - 6 


3.2474 


1.233053-10" 8 


3.5061 


1.059795- 10" 6 


3.4845 


2.343414-10~ 6 



TABLE IV: The same as in Table II but for a different set of nuclei: 154 Nd |2g], 156 Nd |2fl|, 156 Sm 126Q, 158 Sm 
The last energy level of 154 Nd is uncertain and thereby it was not involved in the fitting procedure. 



I60 r 



Gd [20], 162 Gd L2 1 



In Tables II and III are given the results for some isotopes of Th, U, Pu and Cm. Except 
for 228 Th, these isotopes are characterized by large values for the deformation parameter 
d. As we already mentioned, 228 Th has features which are specific to the triaxial nuclei. 
We notice the small values for the r.m.s. obtained in these cases. 

In Tables IV, V and VI are studied the deformed nuclei belonging to the isotopic chains 
of Nd, Sm, Gd, Dy, Er, Yb, Hf, W, Os, respectively. The first six situations can be viewed 
as deformed branches of the nuclear phase transition SU(5) — > SU(3) while the last 
three as the deformed branches of the nuclear phase transition 0(6) — > SU(3). The 
nuclei presented in Table VII, VIII and IX are characterized by small d and moreover 
they satisfy the 0(6) (Table VII), the SU(5) (Table VIII) and X(5) (Table IX) symmetries, 
respectively. For all these nuclei, the defining equation (|3.12|) has been used. One notices 
that this formula for the near vibrational picture describes the excitation energies better 
than the generalized HL formula. This is reflected by the relative r.m.s. values. 
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162 Dy 


164 Dy 


166 Er 


172y b 


174y b 


176 Hf 




Exp. Th(l) Th(2) 


Exp. Th(l) Th(2) 


Exp. Th(l) Th(2) 


Exp. Th(l) Th(2) 


Exp. Th(l) Th(2) 


Exp. Th(l) Th(2) 


2+ 
4+ 
6+ 
8+ 
10+ 

1 9 + 

14+ 
16+ 
18+ 
20+ 
22+ 
24+ 


80.66 80.43 80.59 
265.66 265.25 265.66 
548.52 548.23 548.79 
920.50 921.14 921.53 
1374.80 1374.98 1374.8 


73.39 73.36 73.37 
242.23 242.21 242.23 
501.32 501.34 501.41 
843.68 843.68 843.79 
1261.3 1261.19 1261.3 


80.58 81.00 80.54 
264.99 266.16 264.99 
545.45 547.04 545.56 
911.21 912.25 911.35 
1349.64 1349.13 1349.64 


78.74 78.68 78.71 
260.27 260.18 260.27 
539.98 540.01 540.14 
912.12 912.15 912.27 
1370.07 1370.11 1370.07 


76.47 76.48 76.48 
253.12 253.13 253.12 
526.03 525.93 525.83 
889.93 889.26 888.97 
1336.00 1336.63 1336.0 


88.35 88.08 88.23 
290.18 289.78 290.18 
596.82 596.78 597.32 
997.74 997.97 998.37 
(1481.07) 1481.13 1481.07 


LyKJL.o I7UI.U lyuU-U'* 

2492 2491.26 2489.81 
3138 3138.83 3138.0 
3838 3837 79 3839 85 
4583.15 4591.79 
5370.65 5391.23 
6196.67 6236.35 


1 7A^ Q 1 7A^ 1A 1 7A^ 7£ 

2289.6 2289.65 2289.52 
(2886.0) 2885.99 2885.77 
13528 71 3528 66 3528 7 
(4212.3) 4212.4 4213.49 
4932.68 4936.22 
5685.62 5693.78 


2389.4 2388.37 2390.46 
(2967.4) 2968.66 2967.4 
3576 94 3566 42 
4205.52 4177.02 
4847.84 4789.64 
5498.26 5355.47 


ncifl7 A9.\ 1Qfl7^7 1QP17 79 
\Ly\j/ .*±o) Ly\j/ .OA Ly\J/ -ZZ 

(2518.7) 2518.62 2518.17 
(3198.4) 3198.41 3198.4 
"W42 58 3944 38 
4747.52 4753.51 
5610.15 5623.98 
6527.85 6554.61 


('X £A1 "l 1 R£1 77 1 RAD'}'* 
\LOOi) loDl.Z/ loDU.Za 

(2457) 2456.66 2455.3 
(3117) 3116.76 3115.52 
f 38361 3836 14 3836 
(4610) 4610.05 4612.7 
5434.34 5442.38 
6305.41 6322.51 


LZUo4. 0/ ) ZXiO^.Ly Z.\JOO.OD 

(2646.6) 2646.49 2645.6 
(3308.0) 3308.11 3308.0 
1401 81 401 9 401 3 56 
4747.85 4756.41 
5513.01 5531.95 
6301.26 6336.55 


t + 

E 2 + 

X 


3.29 

0.48 1.10 


3.30 

0.07 0.39 


3.29 

1.14 0.41 


3.31 

0.05 0.22 


3.31 

0.33 1.18 


3.28 

0.23 1.07 


A 1 
A 2 
d 


a 
b 

c 


530.0578 6153.5266 
4.4759 0.0043870 
3.4126 1.171720-10 ~ 6 


607.1406 7551.6774 
3.1089 0.003252137 
3.5835 3.4763788- 10" 7 


850.8099 9490.2037 
0.0142 0.002844946 
3.501835 -6.852260-10" 7 


383.6380 6605.21025 
6.6422 0.00398587 
3.4108 1.605260-HT 6 


587.6927 8756.3360 
4.5844 0.00292061 
3.7067 5.919085- 10 ~ 7 


753.6198 6943.5586 
2.2576 0.00426081 
3.4441 2.818646-KT 7 


ABLE V: The same as in Table II but for a different set of nuclei: 162 Dy [29], 164 Dy |30], 166 Er |31], 172 Yb |32], 174 Yb |33], 176 Hf |3J 





182 |y 


186 


17S Os 


180 Os 


186 0s 




Exp. 


Th(l) 


Th(2) 


Exp. 


Th(l) 


Th(2) 


Exp. 


Th(l) 


Th(2) 


Exp. 


Th(l) 


Th(2) 


Exp. 


Th(l) 


Th(2) 


2" 




100.11 


99.55 


99.82 


122.63 


119.61 


121.30 


131.6 


125.48 


129.07 


132.11 


130.643 


129.82 


137.16 


133.67 


134.65 


4+ 


329.43 


328.81 


329.43 


396.55 


393.36 


396.55 


397.7 


388.23 


397.7 


408.62 


409.45 


408.62 


434.09 


432.25 


434.09 


6+ 


680.50 


680.92 


681.50 


809.25 


809.50 


811.02 


761.00 


751.75 


762.39 


795.08 


795.92 


796.33 


868.94 


870.76 


871.64 


8+ 


1144.4 


1146.29 


1146.07 


1349.20 


1352.17 


1349.20 


1193.80 


1188.34 


1194.95 


1257.44 


1256.25 


1257.44 


1420.94 


1422.74 


1420.94 


10+ 


1711.90 


1713.7 


1711.9 


2002.4 


2003.71 


1998.63 


1681.6 


1685.41 


1681.6 


1767.57 


1766.71 


1766.95 


2067.95 


2065.64 


2061.88 


12 


+■ 


2372.3 


2371.37 


2367.76 


(2750.9) 


2746.52 


2750.9 


2219.4 


2230.28 


2216.81 


2308.71 


2310.46 


2308.71 


2781.26 


2781.84 


2781.26 


14+ 


(3112.3) 


3107.73 


3103.24 


(3562.4) 


3564.28 


3601.12 


2804.3 


2812.15 


2799.12 


2875.0 


2874.35 


2872.48 


3557.7 


3557.76 


3571.28 


16+ 


(3909.2) 


3912.02 


3909.2 




4442.48 


4546.91 


(3429) 


3419.0 


3429.0 




3446.79 


3451.67 




4382.72 


4427.8 


18+ 


(4747.1)? 


4774.46 


4777.95 




5368.52 


5587.55 


4035.32 


4107.78 




4016.34 


4041.91 




5247.84 


5349.05 


20+ 




5686.41 


5703.17 




6331.66 


6723.24 




4637.75 


4837.1 




4570.48 


4640.25 




6145.34 


6334.65 


22+ 




6640.27 


6679.82 




7322.68 


7954.72 




5183.34 


5618.61 




5094.27 


5244.62 




7067.97 


7385.08 


24 


f 




7629.42 


7703.93 




8333.72 


9282.96 




5554.72 


6453.89 




5568.28 


5853.54 




8008.58 


8501.28 




+ 
T 


3.29 






3.23 






3.02 






3.09 






3.16 






X 




2.16 


10.84 




2.74 


14.73 




8.34 


2.33 




1.15 


1.39 




1.99 


5.79 


A x 


o 


1106.3961 


9719.4346 


1193.1479 


3759.9451 


217.4255 


1297.2467 


398.5784 


2299.2478 


539.1025 


3350.11 


A 2 


6 


1.7794 


0.00343844 


1.0740 


0.0108515 


8.5931 


0.0345814 


5.1661 


0.01935326 


6.7316 


0.0136108 


d 


c 


3.7663 


4.179487- 10 - 7 


3.5059 


1.256224- 10- 5 


2.1840 


3.875729-10- 5 


2.4119 


-1.080998 10- 7 


2.7518 


9.3005-10- 6 



TABLE VI: The same as for Table II but for a different set of nuclei: 182 W 123, 186 W 136Q, 178 Os 1370, 180 Os [38], 186 Os |3ffl. The last ener^ 
level of 182 W is uncertain and therefore it was not considered in the fitting procedure. 





170 W 


174 Os 


178 Os 


176 Pi 


178p t 


180 Pt 




Exp. 


Th(l) 


Th(2) 


Exp. 


Th(l) 


Th(2) 


Exp. 


Th(l) 


Th(2) 


Exp. 


Th(l) 


Th(2) 


Exp. 


Th(l) 


Th(2) 


Exp. 


Th(l) 


Th(2) 


2" 




1 ^£ 77 
LDO./Z. 


1 A 7 Q1 

14/ .71 


101. 4z 


lOo.OU 


1 ^9 09 


1 AQ 9^ 


1 'SI £ 


11Q 97 

117.Z/ 


1 9Q 07 

1Z7.U/ 




Zoa.ol 


970 Q1 
ZZU. "1 


/1 7n 1 \ 
(1/U.l) 


1 AQ 9Q 
14:7. Z.y 


1 AQ 7^ 




1 AX. (Sk 


1 XSK AO 


4+ 


462.33 


472.55 


462.33 


435.00 


432.32 


435.0 


397.7 


404.65 


397.7 


564.1 


549.32 


564.1 


(427.1) 


426.55 


427.1 


410.74 


414.33 


410.74 


6+ 


875.53 


881.93 


877.80 


777.63 


778.99 


796.32 


761.00 


769.78 


762.39 


905.6 


918.39 


952.62 


(764.6) 


779.04 


778.67 


757.07 


765.22 


768.11 


8- 




1363.40 


1356.16 


1363.4 


1171.93 


1177.95 


1206.54 


1193.8 


1196.31 


1194.95 


1305.7 


1332.77 


1372.98 


(1177.6) 


1194.67 


1189.99 


1181.50 


1185.52 


1190.82 


10 


+ 


1901.5 


1888.66 


1902.46 


1617.5 


1624.47 


1656.71 


1681.6 


1678.19 


1681.6 


1764.8 


1789.52 


1824.74 


(1660.4) 


1669.39 


1660.4 


1674.28 


1671.64 


1674.28 


12+ 


2464.3 


2476.63 


2488.05 


2113.8 


2116.53 


2144.95 


2219.4 


2212.86 


2216.81 


2277.0 


2287.36 


2310.5 


(2207.6) 


2201.46 


2192.86 


2229.2 


2222.0 


2219.37 


14+ 


(3118.0) 


3118.64 


3118.0 


2656.3 


2653.1 


2672.04 


2804.3 


2799.02 


2799.12 


2833.5 


2825.61 


2833.5 


(2811.9) 


2789.99 


2790.85 


2841.5 


2835.83 


2828.63 


16+ 


(3815.9) 


3813.9 


3792.44 


3239.8 


3233.59 


3239.8 


(3429) 


3435.95 


3429.0 


3423.8 


3403.88 


3396.93 


(3457.5) 


3434.47 


3457.5 


3504.8 


3512.67 


3504.8 


18+ 




4561.93 


4512.55 


3861.8 


3857.66 


3850.29 




4123.2 


4107.78 


4041.80 


4021.95 


4003.64 


(4107.9) 


4134.61 


4195.33 


4252.8 


4252.27 


4250.34 


20+ 




5362.43 


5279.94 


4524.9 


4525.07 


4505.56 




4860.49 


4837.1 


4690.40 


4679.66 


4656.13 




4890.21 


5006.32 




5054.43 


5067.32 


22+ 




6215.19 


6096.38 


5233.0 


5235.69 


5207.47 




5647.65 


5618.61 


5377.0 


5376.91 


5356.52 




5701.14 


5892.01 




5919.06 


5957.41 


24 


+ 




7120.05 


6963.58 


5987.10 


5989.4 


5957.69 




6484.53 


6453.89 


6106.60 


6113.64 


6106.6 




6567.32 


6853.56 




6846.07 


6921.93 


26 


+ 




8076.93 


7883.16 


6786.1 


6786.12 


6757.65 




7371.04 


7344.39 


6878.6 


6889.78 


6907.87 




7488.66 


7891.85 




7835.39 


7961.93 


28 


+ 








7628.4 


7625.81 


7608.61 


























30+ 








8511.6 


8508.41 


8511.60 


























32+ 








9429.7 


9433.9 


9467.53 


























K 

El 


+ 


2.95 






2.74 






3.02 






2.14 






2.51 






2.68 






X 




8.64 


11.98 




4.01 


23.78 




7.20 


2.33 




17.29 


37.15 




17.56 


31.75 




6.15 


8.80 


Ai 


a 


301.1197 


1404.2605 


253.6741 


801.1928 


260.4042 


1297.25 


277.2380 


407.296 


222.5347 


551.9787 


198.9086 


761.038 


A 2 


b 


6.4525 


0.03768205 


5.3413 


0.0673811 


6.1596 


0.0346 


4.9031 


0.2278488 


6.8660 


0.10072468 


7.7625 


0.0650988 


d 


c 


1.6390 


3.305334-10" 5 


1.5935 


8.260490- 10" 5 


1.6566 


3.88- 


10" 5 


1.4494 


3.279647- 10 " 4 


1.5878 


3.293427- 10" 4 


1.5869 


1.71825-10" 4 



TABLE VII: The same as in Table II but for a different set of nuclei: J7n W [39], l7l Os [33], 7s Os [37], J7, 'Pt [34], J7N Pf: [37j], 180 Pt S Also 



the predictions Th(l) are obtained with the expression (|3.12|) corresponding to the expansion characterized by d- small. 





108 Te 


150 Sm 


152 Gd 


154 Dy 




Exp. 


TU M \ 

Th(l) 


TWA 

Th(2) 


Exp. 


TU/1 \ 

Th(l) 


Th(2) 


Exp. 


Th(l) 


TU /OX 

Th(2) 


Exp. 


TU/1 \ 

Th(l) 


TU ZO\ 

Th(2) 


2+ 


625.20 


605.01 


559.90 


333.86 


305.14 


300.10 


344.28 


336.40 


312.00 


334.58 


337.10 


299.71 


4~ 




1289.00 


1296.67 


1289.0 


773.24 


774.07 


773.238 


755.40 


759.65 


755.40 


747.04 


756.47 


747.04 


6- 




2047.9 


2075.43 


2083.78 


1278.75 


1300.12 


1301.78 


1227.38 


1233.74 


1239.98 


1224.08 


1227.93 


1241.06 


8+ 


2945.0 


2936.12 


2945.0 


1836.87 


1857.81 


10C7 ,4 I 

1857.43 


1746.78 


1747.55 


1754.67 


1747.82 


1741.62 


1764.45 


10 


+ 


3886.2 


3876.53 


3883.14 


2433.00 


2438.50 


2433.0 


2300.4 


2297.09 


2300.4 


2304.3 


2293.98 


2315.38 


12+ 


4909.10 


4895.60 


4909.10 


3048.20 


3038.44 


3027.07 


2883.80 


2880.59 


2880.62 


2892.60 


2883.42 


2895.80 


14 


+ 


5980.3 


5992.78 


6032.53 


(3675.70) 


3655.73 


3640.20 


3499.20 


3497.12 


3499.20 


3508.60 


3509.09 


3508.60 


16 


+ 




7167.72 


7261.50 


4305.90 


4289.29 


4273.71 


4142.70 


4146.16 


4159.89 


4172.70 


4170.53 


4156.89 


18 


+ 




8420.22 


8602.58 


4929.20 


4938.46 


4929.20 




4827.39 


4866.08 


4868.60 


4867.44 


4843.65 


20 


+ 




9750.14 


10061.10 


(5592.8) 


5602.83 


5608.38 




5540.59 


5620.76 


5589.90 


5599.64 


5571.61 


22+ 




11157.40 


11641.10 




6282.11 


6312.92 




6285.64 


6426.51 


6349.9 


6367.01 


6343.24 


24 


+ 




12641.90 


13346.10 




6976.11 


7044.47 




7062.43 


7285.57 


7160.70 


7169.45 


7160.70 


26+ 




14203.70 


15178.60 




7684.69 


7804.56 




7870.89 


8199.80 


8027.50 


8006.90 


8025.91 


El 


+ 


2.06 






2.32 






2.19 






2.23 






X 




15.74 


34.41 




16.45 


22.54 




4.45 


14.00 




9.74 


15.61 


Ai 


a 


580.6812 


460.8427 


532.8764 


635.7440 


409.7560 


398.91 


395.7184 


481.5077 


A 2 


b 


9.6296 


0.6424835 


1.7549 


0.194008 


3.9253 


0.360103 


4.3456 


0.270578 


d 


c 


1.2510 


1.419787-10" 3 


1.5563 


7.927202- 10" 5 


1.4306 


4.278524- 10" 4 


1.4135 


2.459865- 10" 4 



TABLE VIII: The same as for Table VII but for a different set of nuclei: 108 Te j4(]]/ 50 Sm |4l|], 152 Gd |42(], 154 Dy 



In Table X we present the results for 150 Nd, 152 Sm and 154 Gd. These nuclei play the 
role of critical points for the phase transitions SU(5) — > SU (3) in the respective isotopic 
chain |13-|l6]. The first two have been presented also in Table IX where they have been 
described by a formula obtained for a near vibrational regime. However, one expects 
that for nuclei close to the critical point the other formula using an asymptotic expansion 
in terms of 1/x works as well. This is actually confirmed by the data presented in Table 
X for the first two nuclei. The isotope 154 Gd is supposed to satisfy the so called X(5) 
symmetry [13]. Our results show that the ground band energies of this critical nucleus 
is described quite well by the compact formulas (|3.9|) , (|3.12[) . 
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150 Nd 


152 5m 


156 Dy 


,r 


Exp. 


Th(l) 


Th(2) 


Exp. 


Th(l) 


Th(2) 


Exp. 


Th(l) 


Th(2) 


2' 


V 


130.21 


122.59 


126.16 


121.78 


107.02 


118.00 


137.77 


83.82 


128.88 


4- 


V 


381.45 


386.14 


381.45 


366.48 


372.03 


366.48 


404.19 


398.35 


404.19 


Q- 


V 


720.4 


726.86 


722.89 


706.88 


719.04 


709.74 


770.44 


797.667 


786.22 


8~ 


V 


1129.7 


1130.85 


1129.70 


1125.35 


1131.97 


1125.35 


1215.61 


1252.12 


1242.82 


10 


+ 


1599.00 


1593.52 


1595.26 


1609.23 


1605.59 


1603.21 


1725.02 


1752.12 


1753.75 


12 


+ 


(2119.00) 


2112.90 


2119.0 


2148.51 


2137.63 


2140.1 


2285.88 


2293.59 


2307.61 


14 


+ 


(2682.50) 


2687.99 


2702.57 


(2736.01) 


2726.97 


2736.01 


2887.82 


2874.50 


2898.42 


16 


+ 




3318.23 


3348.22 


(3362.0) 


3372.98 


3392.27 


3523.3 


3493.72 


3523.30 


18 


+ 




4003.30 


4058.16 




4075.29 


4110.63 


4178.10 


4150.56 


4181.18 


20 


+ 




4742.97 


4834 33 




4833 64 


4892.87 


4859 00 


4844.58 


4871.98 


22 


+ 




5537.10 


5678.39 




5647.87 


5740.62 


5573.00 


5575.48 


5596.23 


24 


+ 




6385.58 


6591.67 




6517.87 


6655.29 


6328.70 


6343.07 


6354.73 


26 


+ 




7288.35 


7575.26 




7443.56 


7638.09 


7130.30 


7147.20 


7148.47 


28 


+ 




8245.35 


8630.02 




8424.87 


8690.01 


7978.50 


7987.76 


7978.50 


30+ 




9256.53 


9756.66 




9461.76 


9811.89 


8875.90 


8864.66 


8845.86 


K + 
E l + 


2.93 






3.01 






2.93 






X 




5.61 


7.92 




9.84 


11.43 




23.86 


18.45 


A x 


a 


215.8125 


2867.41 


221.5088 


1187.15 


348.6711 


1913.4867 


A 2 


b 


6.7513 


0.0513 


6.9223 


0.0344 


4.4983 


0.0231439 


d 


c 


1.6321 


1.17-10- 4 


1.6640 


6.IMO- 5 


1.7089 


1.051357-10- 5 



ABLE IX: The same as in Table VII but for a different set of nuclei: i&u Nd [41], i&2 Sm [42], 
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150 Nd 


152 5m 


154 Gd 


,r 


Exp. 


Th(l) Th(2) 


Exp. 


Th(l) Th(2) 


Exp. 


Th(l) 


Th(2) 


2+ 


130.21 


123.45 126.16 


121.78 


119.14 118.00 


123.07 


108.23 


117.32 


4+ 


381.45 


374.67 381.45 


366.48 


364.57 366.48 


370.99 


346.04 


370.99 


6+ 


720.40 


717.10 722.89 


706.88 


704.82 709.74 


717.65 


689.42 


728.47 


8+ 


1129.70 


1130.89 1129.70 


1125.35 


1123.09 1125.35 


1144.43 


1117.23 


1162.21 


10+ 


1599.00 


1603.40 1595.26 


1609.23 


1608.64 1603.21 


1637.04 


1614.03 


1654.28 


12+ 


(2119.00) 


2123.44 2119.00 


2148.51 


2151.71 2140.10 


2184.67 


2168.74 


2194.41 


14+ 


(2682.5) 


2677.98 2702.57 


(2736.01) 


2740.64 2736.01 


2777.30 


2772.90 


2777.30 


16+ 




3248 1 2 3348 22 


(3362 0) 


3357 92 3392 27 


3404.44 


3419.34 


3400 53 


18+ 




3798.63 4058.16 




3968.95 4110.63 


4087.10 


4101.33 


4063.38 


20+ 




4221.54 4834.33 




4444.16 4892.87 


4782.30 


4811.93 


4766.07 


22+ 




5678.39 




5740.62 


5519.50 


5543.34 


5509.31 


24+ 




6591.67 




6655.29 


6294.10 


6286.21 


6294.10 


26+ 




7575.26 




7638.09 


7055.50 


7028.36 


7121.53 


< 


2.93 




3.01 




3.01 






X 




4.83 7.92 




2.93 11.43 




21.22 


21.77 


Ai a 


151.3419 


867.4073 


126.4776 


1187.1506 


274.9237 


1950.92 


A 2 b 


9.5423 


0.0513027 


10.3917 


0.034413709 


7.4351 


0.02057356 


d c 


2.0284 


1.171278-10- 4 


2.0136 


6.109010-10- 5 


2.4767 


1.254659- 1Q- 5 



ABLE X: The same as in Table II but for a different set of nuclei: 150 Nd 141D, 152 Sm ||42 
123- 
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In the last table (Table XI) we present two nuclei which satisfy the symmetry E(5). 
These are described with the close formulas (|3.12|) and (|2.54|) . We remark that also in this 
case the r.m.s. values are small. Experimental data for these nuclei were considered up 
to the angular momentum were the first backbending is showing up. 





104 Ru 


W2 Pd 


J 


IT 


Exp. 


lh(l) 


lh{2) 


Exp. 


lh(l) 


lh(2) 


2- 


h 


358.03 


348.13 


358.03 


556.43 


563.02 


520.03 


4" 


h 


888 49 


901 59 


910 87 


1275.87 


1280.12 


1275.87 


6" 


h 




1 ££1 OQ 
lOOl.ZV 




01 1 1 3£ 


01 nn 1 
zlUU.lZ 


011/1 71 
Z114./ 1 


8' 






2303.22 


2287.30 


301 3 06 


3006.50 


3019.06 


10 


+ 


3111 80 


3119.23 


3111.80 




3993.31 


3992.71 


12 


+ 




4005.75 


4039.35 


^0^ 10 


5057.87 


5043.55 


14+ 




4960.97 


5078.09 


6179.80 


6198.80 


6179.80 


16 


+ 




5983.88 


6234.53 


7428.80 


7415.32 


7408.97 


18 


+ 




7073.86 


7513.64 




8706.93 


8737.58 


20+ 




8230.5 


8919.18 




10073.30 


10171.10 


*f 
E t 


2.48 






2.29 






X 




11.33 


17.84 




9.88 


15.41 


A x 


a 


522.4864 


562.6938 


649.8648 


710.6779 


A 2 


b 


8.1960 


0.2738559 


9.2391 


0.329923 


d 


c 


1.5466 


9.51995-10- 4 


1.4174 


5.317905-10- 4 



TABLE XI: The same as in Table VII but for a different set of nuclei: 104 Ru 
nuclei obey the E(5) symmetry. 



102t 



Pd ||44Q. These 
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V. CONCLUSIONS 



In this Section we summarize the main results obtained in this paper. By a de- 
quantization procedure we associated to a quantum mechanical Hamiltonian which is 
quadratic in the quadrupole bosons, a time dependent classical equation. The classi- 
cal Hamiltonian has a separated form, i.e. is a sum of a kinetic and a potential energy 
terms. The later one is not depending on momenta and is of Davidson type. We may say 
that our procedure proves the classical origin of the Davidson potential. The centrifugal 
term is determined by a pseudo-angular momentum associated to the intrinsic coordi- 
nates. It is worth mentioning that the constraint for the angular momentum in the lab- 
oratory frame yields a differential equation which is connected to that one correspond- 
ing to the energy conservation which results in obtaining a specific angular momentum 
dependence for the quantal energy. Actually the expression obtained generalizes the 
Holmberg-Lipas formula, involving under the square root symbol a J 2 ( J + l) 2 term. 

A similar expression was obtained by one author (A. A. R.) within the coherent state 
model (CSM) for a large deformation regime. Another compact expression was pro- 
posed by CSM for the near vibrational regime, i.e. small nuclear deformation. One of 
the targets of this paper was to prove that the two compact expressions provided by 
CSM are able to describe the ground state energies for deformed, near vibrational and 
transitional nuclei. By matching the two expressions one obtains a unitary description 
for nuclei satisfying different symmetries or, with other words, belonging to various nu- 
clear phases. Similar goal is touched by using a square root formula, a generalization of 
the Holmberg-Lipas formula, obtained on the base of a semiclassical description. 

These descriptions are used for a large number of nuclei (44). The agreement between 
results and experimental excitation energies is very impressive. The agreement quality 
is judged by the small r.m.s. values of discrepancies. 

As a final conclusion we may say that the CSM procedure is able to describe in a real- 
istic fashion the ground state energies for nuclei of different nuclear phases. An alterna- 
tive description is given by a square root formula derived as approximate eigenvalue of 
a quadratic Hamiltonian in quadrupole bosons subject to a constraint due to the angular 
momentum conservation. 

Acknowledgment. This work was supported by the Romanian Ministry for Educa- 
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VI. APPENDIX A 



Using the equations of motion for the conjugate variables, one can prove that 



£ 3 = o, m = 0, 

where £3 is defined by the following expression: 

h, , h 2 



(A.l) 



(A.2) 



£3 = ^(qiP2 - qzPi) = -j;r 9, 

and has the significance of the third component of the angular momentum defined in 
the phase space, spanned by the coordinates (gi,pi, q2,P2)- The other two components 
are: 

h h 
A = j{(Qi +Pi -?2 -Pl)> ^2 = -^{qm + PiPi)- (A3) 

Indeed, one easily check that 

{Ci,C k } = he ikj Ci, (A A) 

where {, } denotes the Poisson bracket while e ik j the antisymmetric unit tensor. In virtue 
of Eq. (|A.4|) the set of functions C k with the Poisson brackets as multiplication operation, 
form a classical SU C (2) algebra. Moreover, they could be obtained by averaging with 
I \&), the generators L k 

A = (* I = 1,2,3, 

of a boson SU b (2) algebra defined with the boson operators &£, b\_ 2 , as: 

h 



(A.5) 



L 

L 3 



h 



2V2 
h 

2^2i 



2blb -(bl + bl 2 )(b 2 + b^ 
bl(b 2 + 6_ 2 ) + (4 + bl % )k 



bl(b 2 + b^)-(bl + bl 2 )b G 



(A.6) 



The equation (2.15) and the correspondence between commutators and Poisson brack- 
ets [, ] — > 4{, }, define a homeomorphism of the boson and classical algebras generated 
by {I/fc}fc=i,2,3 and {£fc}fc=i,2,3 respectively. Note that the boson SUb{2) algebra does not 
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describe the rotations in the real configuration space but in a fictitious space. The conser- 
vation law expressed by (|A.1|) is determined by the invariance against rotation around 



the 3-rd axis in the fictitious space mentioned above: [H, L 3 ] = . Since the classical 
system is characterized by two degrees of freedom and, on the other hand, there are two 
constants of motion 

H = E , £ 3 = L, (A.7) 
the equations of motion are exactly solvable. 

VII. APPENDIX B 

By direct calculations we can check that the overlap integral if and its first and 
second derivatives satisfy the following differential equation: 

<Plf X -Zdlf 2x 2 + J(J+l) (0) 
dx 2 2x dx 4x 2 J _U ' (OA) 



By a suitable change of function this equation can be brought to the differential equation 



characterizing the hypergeometric function of the first rank. Thus the final result for if 



is: 

/(°) = ( J!)2 (M 2 )ie~ d ^r l F l ( -(J + 1), J + -; -d 2 ) . (B.2) 
J (|)!(2J+ l)! v ; V2 2' 2 J v 7 

This expression is further used for describing both the asymptotic and vibrational be- 
havior for the excitation energies in the ground band. Indeed, in the asymptotic region 
of d, the hypergeometric function behaves like: 

lFl(a ' c; z) = Y§) eZza ' c[l + ^M" 1 )]- (B3) 

Due to this expression the dominant term of I^ 1 is: 



!f ~ i-. (B.4) 



This expression suggests for I^p , in the asymptotic region, the following form: 

if = e x J2 A n x~ n - (B.5) 



n=l 
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Inserting this expression into the above differential equation one obtains the recursion 
relation for the expansion coefficients A k : 



A n+1 = ^(2n + J)(2n-J-l). 
bn 



(B.6) 



The leading term (|B.4[) gives A% = | and then (|B.6|) determines the whole set of the 
expansion coefficients. In this way we obtain for the ratio cPlj/Ij the expression: 



x- 



'I 

r(0) 



x — 1 — 



1 
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3x 9x 2 27x 3 



1 5 

+ 



13 



6a; 18x 2 18x s 



J(J + 1) 



54x 3 



j 2 (J+l) 2 + 0(x- 



(B.7) 



The convergence in terms of x for the excitation energy may be improved in two steps. 
First we write the differential equation for if" 1 in a different form: 



r(l)' 



r(i)' 



X X 



+ x- 



x — 1 



X- 



r(o) / 1 I r(o) / 2 \ r(°) 



2x 2 + + 

4 



0. 



(B.8) 



The derivative (^707 J is further calculated by using (B. 7)(B. 7)(B. 7)(B. 7)(B. 7)(B. 7)(B. 

7 (i) 

7)) and thus the above equation becomes a second degree algebraic equation for x^j. 



Solving this equation one obtains: 



x- 



r(l) 

j 

r(0) 



x - 2 



(B.9) 



where we used the notation: 



G J= M*-2) + (J + ±) -^(3 + f + £) 

+i(i + g + S)^+i)-S(^+i) 2 - 



(B.10) 



Concerning the near vibrational regime the final expression for energies is obtained in 
two steps. First one derives the vibrational limit of the k-th derivative: 

r( i)\ W 



lim I d 



I 1 J 



r(0) 



(2J + 3)' 



(J + l)(J + 2) 
+9 2J + 5 (fe + 9 



J 
2 

<5fc,3 



2J + 7 



, = 0, 1,2,3. 



(B.ll) 
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j(l) 

Then, truncating the Taylor expansion of x-fa, around the point x = 0, at the third order 
one obtains: 

if J J 9 (J+l)(J + 2) 2 

^ ~ 2 + 2(2J + 3) ;r+ 2 2J + 3) 2 (2J + 5) ;r 

+ 2 (2J + 3)3(2J + 5)(2J + 7) +Ul j - { ' 
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